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Summary 

A mathematical model for a general multibody flexible spacecraft is obtained. The generic 
spacecraft considered consists of a flexible central body to which a number of flexible multibody 
structures are attached. The coordinate systems used in the derivation allow effective decoupling 
of the translational motion of the entire spacecraft from its rotational motion about its center of 
mass. The derivation assumes that the deformations in the bodies are only due to elastic motions. 
The dynamic model derived is a closed-form vector-matrix differential equation. The model 
developed can be used for analysis and simulation of many realistic spacecraft configurations. 

1 Introduction 

A class of next-generation spacecraft is expected to inlude nonlinear, multibody, flexible space 
systems. Some of the current spacecraft can also be catagorized under this class. Examples 
of such systems are [1] : satellites with flexible appendages, such as solar arrays and antennas, 
space-shuttle with remote manipulator system (RMS), and flexible space platforms with multiple 
articulated payloads. Mathematical modeling of these systems is quite complex. This problem 
has been addressed in the existing literature (e.g., see [2]); however, in this paper, a different 
approach is taken to derive the equations of motion which yields compact closed-form equations of 
motion. The derivation uses modeling techniques similar to those used in the robotics literature 

(e.g., see [3], [4]). The formulation is relatively general and can be used for a large class of 
spacecraft. 

First, for the sake of completeness, some of the mathematical aspects of modeling rotat- 
ing systems are summarized in the section on mathematical preliminaries. Next, the kinematic 
equations, i.e., the position, velocity, and acceleration relations, for a particle mass of the sys 
tern aie obtained. Once the kinematic equations are derived, the dynamics of the system can 
be modeled by using various methods; for example , the Newtonian approach, the calculus of 
variations approach, and the Lagrangian approach. The equations of motion derived using any 
of these methods are equivalant; however, the Lagrange- Euler formulation is used here since it is 
an energy-based approach (i.e, it uses scalar formulation) and is easy to work with compared to 
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Newtonian approach which deals with vector quantities. It is assumed in the derivation that the 
bodies deform only due to the elastic motion. However, any other deformations such as thermal 
deformations can be easily included in the formulation with some modifications in the potential 
energy function. 

In deriving the kinematic equations for the chain of multiple flexible bodies, the coordinate 
systems become an important element of the derivation. A large part of kinematics deals with 
the coordinate transformations used to represent the position and orientation of the body. In 
view of this, we will begin with a study of the operations of translation and rotation, and the 
transformations which are used to represent these motions. 


2 Mathematical Preliminaries 

2.1 Rotations 

C onsider a rigid body as shown in Figure 1, to which a body-axis system, i.e., a body-fixed 
coordinate frame, ( X b , Y b , Z b ), is attached. Let ( X u , Y 0 , Z Q ) represent some fixed reference frame 
whose origin is concentric with the body axes system. Our aim is to relate the coordinates of 
a point P on the body in the ( X b ,Y b ,Z b ) frame to those in the ( X 0 ,Y 0 ,Z 0 ) frame. Let i b ,j b ,k b 
denote orthonormal basis vectors in the body frame and i 0 -,joik 0 denote the orthonormal basis 
vectors in the fixed frame. Then, the position vector of point P, f p , in the body frame is given as 

= r bxH + r by j b + r bz k b (1) 

and in the fixed frame is given by 

T'ox^'o ~f* T'oyjo “b T'ozfco (2) 

Since r p and r p are representations of the same vector r p , the relation between the components 
of r p in two different systems can be obtained as follows. 

r bx = h * r p = *6 * (3) 


Similarly, 



and 


n Z = h ■ r b p = k b -r° p 


(5) 


Equations (3)-(5) can be rewritten in a compact form 


as 


“ 


- 

rbx 


r Q x 

V by 

= R b 0 

r oy 

r bz 

- 


r oz 


or, notationally, 


where 


Similarly, the coordinates of 
in the body frame as 


where 





= R b f° 



(6) 


H 

1 0 

* Jo 


K 


R b 0 = 

Jb 

1 0 

Jb * Jo 

Jb 

k Q 

(7) 


h 

2 o 

h ■ jo 

kb 

k 0 


in the 

fixed frame can 

be 

expressed in terms of the coordinates 



n-- 

~ n b r p 



(8) 


1*0 


^ o * Jb 


h 


Rl = 

jo 

• h 

jo * jb 

jo • 

h 

(9) 


Jc 0 

^ b 

k o ‘ Jb 

k Q 

h 



Since dot products are commutative, from equations (7) and (9), we can see that 

R° b = Rf = R b ~ l 

and 

R b o = Rf = Rr 1 


Then, 


R b 0 R° b = RfRl = R bT R b = /, 


( 10 ) 


( 11 ) 


(12) 


where I 3 is the 3 x 3 identity matrix. The transformations, R° b and R b 0 , are called ‘orthogonal’ 
ti ansformations. These transformations are also referred to as ‘orthonormal’ transformations 
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since all column vectors in the transformation matrix are unit vectors in addition to being 
orthogonal. Thus, the transformation matrices can be used to relate the representations of the 
same vector in two different coordinate reference frames. Note that the transformations do not 
change the vector itself but only its representation. Another important thing to be noted is that 
the rows of R% are the direction cosines of i Q , j a , and k a , respectively, in the body coordinate 
frame. 

Matrices R b and can also be interpreted as rotation matrices wherein {ib^jb^b) are or- 
thonormal basis vectors in the final direction of the {i 0 ^j 0 ^o) axis system after rotations about 
the selected coordinate axes. The properties of rotation matrices are developed in the next 
section. 

2.2 Basic Rotation Matrices 

When the rotation matrix represents a change in orientation about any one of the principal 
coordinates of the reference frame, X, V, or Z, it is called as the ‘basic’ rotation matrix. So, if 
a new coordinate system, say B{u,v,w), is obtained by rotation a of the old coordinate system 
, say A(u,v,w), about X axis, then the basic rotation matrix associated with this rotation is 
given by, Rg = R Xt0t . Similarly, the basic rotation matrices associated with rotations, /? about Y 
axis and 7 about Z axis, are given by Ry,p and Rz } 7 , respectively. Referring to Figures (2a)-(2c), 
the basic rotation matrices can be written as 

1 0 0 

Rx,ct = 0 cosa —sina ( 13 ) 

0 sina cosa 

cos /3 0 sin /3 

Ry 4 s =010 ( 14 ) 

—sinj 3 0 cosf 3 

cos 7 —sin*f 0 

Rz,i = sin~j C057 0 

0 0 1 
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The reason these matrices are called basic rotation matrices is because any finite arbitrary ro- 
tation can be achieved by a composition of these matrices. However, since the finite rotations 
are not commutative, the order of multiplication of these matrices during composition is very 
important. 

2.3 Composite Rotations 

As stated previously, any arbitrary finite rotation can be achieved by a composition of the basic 
rotation matrices, i.e., by following a sequence of basic rotations. In obtaining the composite 
lotation matrix there are three different possibilities. The successive rotations can take place 
either about the prinicipal axes of the fixed reference frame, or it can take place about the 
piincipal axes of the rotating frame itself, or a combination of both. The following procedure 
can be followed to obtain a composite rotation matrix. 

When the rotation occurs about any principal axis of a fixed reference frame, premultiply the 
last resultant rotation matrix by the corresponding basic roation matrix and, when the rotation 
occurs about any principal axis of rotating reference frame itself, postmultiply the last resultant 
rotation matrix by corresponding basic roation matrix. 

Let us suppose that the two axes systems, OXYZ (fixed) and oxyz (rotating), are initially 
coincident. Then the rotation matrix will just be an identity matrix, 1 3. Now suppose that oxyz 
undergoes the following sequence of rotations. It rotates about OX axis through an angle ot and 
then rotates through an angle (f> about ox axis. The composite rotation matrix in this case will 
then be given by 

R = Rx,cxhRx,<t> — Rx,aRx,<f> = Rx } (a+<t>) ~ Rx,(a+<f>) (16) 

1 0 0 1 0 0 1 0 0 

0 cos a - sina 0 coscf> -sin<j> = 0 cos(a + <f>) -sin(a + <j>) (17) 

0 sinct cosat 0 sincp cos<f) 0 sin(ot T <f>) cos{a -f <f>) 

2.4 Rotation About an Arbitrary Axis 

In many cases, the rotations of the body-fixed coordinate system take place about an axis other 
than the principal axes of the fixed frame or the rotating frame. In the case when the rotation 
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takes place about an arbitrary axis, the rotation matrix can be obtained as follows. Referring to 
Figure 3, let v = {v x ,v y , v z } T be the unit vector along the axis of rotation and ip be the angle 
through which the rotation takes place. Now obtain the rotation matrix for following sequence 
of operations. Rotate v about OX by angle a which will bring vector v in the XZ plane. Then 
rotating about 0} axis by — /? will align v vector with OZ axis. Then rotate about OZ or v by 
angle ip and then reverse the rotation order to bring vector v back to its original position. This 
sequence of rotation will lead to the following composition of basic rotation matrices. 


Noting that, 


R — Rx^aRy^Rz^Ry-pRx^ 


(18) 


S 2 71 O' 


cos a — 


\l V l + v l \J V l + V, 

P = V x cos/3 = yjv* + V 


sin 


the rotation matrix can be rewritten as [3] 

l ’J(l — cosip) + cosip v x v y {\ — cosip) — v z sinip v x v z {l — cosip) + v y sinip 

1 — cosip) -)- v z sinip v yP — cosip) -f- cosip VyV z {\ — cosip) — v x sinip 

[ tvMl “ cosil>) - Vysinip t-’ y n 2 (l - cosip) + v x sinip v 2 z {\ - cosip) + cosip 


R 


(19) 

( 20 ) 


( 21 ) 


In mathematical notation, these 3x3 rotational matrices are said to belong to SO( 3) space. 
The notation SO{ 3) stands for Special Orthogonal group of order 3 [4]. 


2.5 Properties of Rotational Matrices 

The rotational matrices have some special properties which play an important role in the math- 
ematical modeling of the system. These properties are listed below. 

1) As shown previously, R r = i? _1 , i.e., RR T = / 3 . 

2) The columns of R represent the unit vectors along the principal axes of the rotated coordinate 
frame with respect to the reference frame unit vectors. 

3) Since R R the rows of rotation matrix represent the unit vectors along principal axes of 
reference frame with respect to rotating frame. 
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4) Any row (column) of rotation matrix is orthonormal to any other row (column). This is the 
direct consequence of properties 1 and 2. 

5) If a, b E 7Z 3 , where 7Z 3 is 3— dimensional Euclidean space, then R(a xb) = Ra x Rb where 
symbol x denotes vector cross product. (Note that this equality is valid only for orthogonal 
rotational matrices) 

2.6 Skew-Symmetric Matrices and Cross Product Operator 

A skew-symmetric matrix, <5, has the property: s„ = 0 and = —s Jt for i ^ j. Then, an 
immediate consequence of this property is given by: S + S T = 0. These matrices play an 
important role in the computation of vector-matrix operations involving vectors belonging to 7 Z 3 
space and matrices belonging to 7ZO{3) space. To illustrate, consider a vector cross product axb 
which can be written in terms of vector-matrix multiplication as S(a)b, where S(-) is referred as 
the ‘cross product operator matrix’, and is given by 

0 -a z a y 

S(a) = a z 0 -a x (22) 

“™ CLy Ct j 0 

This cross product operator matrix has some important properties which are discussed below. 
For the sake of simplicity, we will use the following simplified notation in the remainder of the 
paper. 

S(-) = (■) (23) 

2.7 Properties of Cross Product Operator Matrix 

1) Using property (5) of rotational matrices and the definition of <S(-), it can be shown [4] that 
for any vector, v € TZ 3 and R € 7ZO(3 ), 

RS(v)R t = S(Rv) (24) 

2) If Rv,e represents the rotation about an axis aligned with vector v by angle 0 then the derivative 
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of R with respect to 6 is given by [4] 


dRv,e 

de 


= S{v)Rs fi 


(25) 


3) As an obvious consequence of the cross product property, we obtain another property: <S(vi)u 2 = 
— <S(u 2 )ni, where Fi and rJ 2 are 3- vectors. 


3 Mathematical Modeling 

The objective of this section is to derive the equations of motion for nonlinear, multibody, flexible, 
spacecraft in the most general framework possible. 

3.1 Modeling Considerations 

The focus configuration of a generic spacecraft under consideration has a branched geometry with 
a relatively large central body. The system under consideration can be schematically represented 
by the configuration shown in Figure 4. It is assumed that all bodies in the system are flexible. 
The deformations in the bodies are assumed to be due to elastic motions only; however, any 
other deformations such as, due to thermal effects, can also be modeled if required. The system 
model under consideration has cluster configuration. It consists of one central body attached 
to various appendage-bodies to form a branch geometry. For the purpose of derivation the 
following notations are used. Let each body be denoted by where, the first subscript indicates 
the branch the body belongs to and the second subscript indicates the body number in that 
particular blanch. Since the number and the locations of various bodies are arbitrary the system 
configuration is highly general. 

3.2 Coordinate Systems 

Referring to Figure 5, Yj , Z/ is the inertial coordinate system; X cm ,Y cm , Z cm is the coordinate 
system with the origin fixed at the center of mass of the entire spacecraft and is aligned with 
the inertial frame, X c , Y c , Z c is the coordinate frame attached to the central body with the origin 
attached to the center of mass, and Xij,Y{j, Z,y represent the local coordinate system attached 
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to the ?j-th body with the origin located at the point of connection between i(j — l)-th body 
and ij - th body. The motion of each local coordinate system origin, O l3 , is defined with respect 
to the previous local coordinate frame. In order to derive dynamical equations, it is necessary 
to obtain expressions for the kinematic quantities, i.e., the position, velocity, and acceleration of 
the spacecraft. 


3.3 Kinematics of a Spacecraft 

Consider a spacecraft with its instantaneous center of mass located at 0 cm (Figure 5). Then, 
the vector, r cm , determines the inertial position of the spacecraft. In order to decouple the 
translational motion of the entire spacecraft from the rotational motion about its center of mass, 
for the reasons that will be apparent later, the orientation of the frame X cm , Y cm) Z cm is assumed 
to be same as A/, T/, Z/, i.e. the rotation matrix, will be an identity matrix. The vector r cm 
can be expressed in terms of orbital elements [2]. Vectors a c and ~p r represent displacement of the 
center of mass of the central body due to rigid and elastic motions, respectively. Vector r c then 
represents the vector sum of these two vectors, i.e., f c = a c + p c . In addition to the translational 
motion of the entire spacecraft, the spacecraft can undergo the rotational motion about its own 
center of mass. This motion can be characterized by the rotational transformation, R c c m , between 
F cm , Z cm and X c , Y c , Z c . Since, X cm , V cm , Z cm and X/, V/, Z/ are aligned, R \ : m describes the 
orientation of the spacecraft with respect to the inertial frame also. R c c m is generally obtained 
by using Euler rotations. Other rigid body degrees of freedom arise from the interconnections 
between different bodies of the spacecraft, each of which can be described by the transformation 
of the type, R^ J between any two consecutive body frames in the chain. Oij represents the 
origin of the ij-th body frame and its position with respect to i(j — l)-th body frame is defined by 
vector Sij. Also, each S{j is the vector sum of and ~p tJ , where a,j and p- represent rigid body 
and elastic displacements of ij-th frame, respectively. First, the equations for position, velocity, 
and acceleration will be obtained for a representative particle mass dm in the ij-th body, i.e., 
the j-th body in the Tth branch. Referring to Figure 6, the position vector of a particle mass 
dm in ij-th body, in the local reference frame (i.e, ij-th frame), is given by 

*i™ = 2?m + Pdm (26) 
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where <T ^ m and ~pj m represent rigid and elastic displacements of mass dm in ij - th frame. The 
position vector of dm with respect to i(j — l)-th frame is then given by 


u 


—- 5 .., piU~ 

dm ij U'dm 


(27) 


Finally, the position vector of dm in the inertial frame of reference is given by 


( k=j - 1 


^ dm — ^ cr 


+ ^cm r c+( 


fc=l 


(28) 


where R{- is given by 


*4 = ■ ■ ■ r?t" 


(29) 


The velocity of the particle mass dm is given by taking the time derivative of equation (28). 


d d ^ 

v dm = bim = r cm + —(R^Sc) + — ( Rik^ik) + 

k=j- 1 ^ Ar=j — 1 

= Vcm + R'cm^c + (Rik^ik) + Y (Rik^ik) 
k = 1 fc=l 

+ + 


(30) 


Now using properties 2 and 3 (Section 2.5) of the cross product operator, noting that R[ m = / 3 , 
and using the notation (23), the time derivatives in equation (30) can be evaluated as follows. 

s c =a c + p c = Lj c a c + — —a,Jo c + Y $ck<jck (31) 

k k 

where, p c ~ J2k ^ckQck, $ck is the mode shape matrix of the central body, and q c k are the 
generalized coordinates. Similarly, 

S,J = dij + Pij = -I- $ijk<iijk = -dijLOij + Y ^ijkQijk (32) 

k k 

where, p tj = (p^kqijk and being the &-th mode shape matrix and generalized coordi- 

nate for ij-th body, respectively. 

r'„s„ = Ar*?, • • • «g i -% + «r«r, • • • tff'% 

+ •••+ r% 

= (Rr r;, ■ ■ ■ r?/-%p c +R r(Rf r - R‘y- , %p ll 
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(33) 


= {RTRU- + *?(%> ■ • • 

+ + < 34 ) 

r’Jl = r'h&L+^L) 

= tfiisA + E****) 

k 

= -MliZ^ + RljEKli* < 35 ) 

k 

Substituting equations (33) thru (35) in equation (30) we get 


^dm — Vcm Q, C UJ C ~f" ^ ^ ^ck^cA: 

k 

— Rr m \ j 

j=i i=i k 

+ • • • K ( r x) ^c + k? [E(^)]^i + 

j J=1 




•O- 1 ) 0 . ^777. 




+ (nr A), • ■ ■ + htWi ■ • • 

+ ■ •■ + «r«ii • • • 

- r;X^» + ■«« E 


t (j-' 1 )r 7 'i 


(36) 


Simplifying and regrouping equation (36) gives 


Vdn 


V C jn 


+ [-# c + j2(R c c m R« • • • *r%o + 




+ [-«S"s« + flT(E(^ J «)) + + 

i=i 




.cm / nUz-l)^ \ 1 — 


]=z 


+ • • • + + ^(7-i)( 4 j_1) ^) + ^-Di^r 1 ^)]^ 


?cm O \lrt. 


- + E *m* + E fi“(E + K E *&««» 

k 3=1 k k 


(37) 
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Equation (37) can be rewritten in a compact form as 


V dm = v l d T s + vZ = Np (38) 

where and v ^ are the contributions to the velocity vector of particle mass dm from the 

translational and rotational-plus-flexible motion of the spacecraft, respectively, and 


—trans 
v dm 


^ cm 


( 39 ) 


- rot 
'dm 


Np 


(40) 


N = [I 3 N) (41) 

? = {Vcm p} T (42) 


In equations (39)-(42), N is 3 X n matrix, n is the total number of rotational-plus-flexible degrees 
of freedom, and p is n x 1 vector of generalized velocities corresponding to these degrees of freedom. 
Matrix N has the following form: 

N [ Nrijid Nfi ex ] (43) 

where, N rigid is 3x3(/m+l) matrix related to the rigid degrees of freedom and N Jlex is 3xs(/m+l) 
matrix related to the flexible degrees of freedom. For the sake of notational simplicity, and 
without loss of generality, it is assumed that each chain structure has m bodies and that each 
body has s flexible degrees of freedom. Then vector p has the form 


P [w c ,0Ju, • • • ^lmi • • ' • • • u>i m , q c i , • • ■ q cs , j 1U) • • • q i j ky • • • qi ms Y (44) 

The matrices N rigid and Nfi ex are given as follows. 


Nrigid — [N c , Nn , ' ■ ' i N tz , • • ■ , N tJ ] 

where 

N c = -S e + ^(^ • • • tfq -%) + (R c c m ^i • 

3 

N« = + Rr(jZ(Rl^ l3 )) + RT(R^-'L)\ 

;=i 


(45) 


(46) 

(47) 
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(48) 


N„ 

N„ 




j=Z 


(49) 

(50) 


and 

• • • • £ *„« < 51 ) 
- it k . 

Having obtained the expression for velocity, the kinetic energy for the particle mass dm is given 

by 

(52) 


1_ T _ 

dTdm = 2 V dm Vdmdrn 


Then the kinetic energy for the entire spacecraft can be obtained by integrating equation (52), 
i.e., 

T = f dT dm = \ j VjnVdmfJ’dfi, (53) 

Jn 2 Jn 

where, p is the mass density, v dm is as given in equation (37), and fi denotes the spatial domain 
of integration. Substituting equation (52) into equation (53), we get 


vT = l [ {N^) T Nffidn 
2 Jn 

= \ [ (N T N)^ndn 

L JQ 


(54) 


Equation (54) can be rewritten in a compact form as 

T = M(p)f 


(55) 


where, M(p) is the mass-inertia matrix of the system and is given by 


M{p) = I ( N T N)fidn (56) 

M(p ) is symmetric and positive definite matrix. 

In practice, for the most spacecraft applications, a control engineer needs to know only 
rotational dynamic model of the spacecraft since the translational motion of the spacecraft, as 
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a whole, is controlled only periodically by separate thruster-jets whenever reboosting of orbit is 
necessary. In such a case, the kinetic energy of the spacecraft is given by 

T '°' = sX’SVmMl 

= \ [ (Npf NptulSl 

Z JQ 

= \ J a P T {N T N)pndn ( 57 ) 

and the inertia matrix is then given by 

M(p) = J^(N T N)pdn (58) 


M (P ) is also symmetric and positive definite matrix. Using (57) and (58), T rot can be rewritten 
in compact form as 


T rot = 


(59) 


3.4 Potential Energy 

The potential energy of the system could be due to many sources, such as elastic displacement, 
thermal deformation, etc. The deformations due to thermal effects is not considered here, how- 
ever, it can be easily included in the formulation if desired. Thus, it is assumed that the potential 
energy has contribution only from the elastic strain energy. Also, it assumed that the materials 
under consideration are isotropic in nature and that they obey Hook’s law. 

For the isotropic materials obeying Hook’s law, the strain energy differential is given as 

8V = J a T Stdil (60) 

which can be rewritten as 

6V = L * <Hi ( 6i ) 

where is the strain energy density and has the form 

® = a xx^xx + Vyytyy + . . . + Cr yz t yz (62) 
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Now, for materials obeying Hook’s law, following equality holds. 


a 1 = Ee 


( 63 ) 


The strain-displacement relation is given by 


e = T>u 


( 64 ) 


where, u is the general displacement vector and T> is the differential operator defined by relations 


C,j “ 2 


u 


,3 + U j,i + E U 'k!> 
k = 1 


{l,j = 1,2,3} 

The vector u can be expressed in terms of modal coordinates as 


Now, from equation (64) 


Substituting in <5V, get 


u = <$>q 


Se = V$8q 


8V 


f <7 T 8cdQ 
Jn 

= [ e T EV^8qdil 

Ju 

= / q T (V<b) T EV$8qdn 

Jo. 

= q T [ (D $) T EVQdttSq 
Ju 

= q T K8q 

where, K is called the stiffness matrix of the system and is given by 

I\ = [ (D $) r EVQdVl 
Jo 

The potential energy of the system is, then given by 


(65) 


(66) 


(67) 


( 68 ) 


(69) 


1 T 

1 / = -q T I<q 


( 70 ) 
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3.5 Equations of Motion 

As previously stated, without the loss of generality, only the rotational motion of the spacecraft 
is considered in deriving the dynamic model of the spacecraft. If necessary, the translational 
motion can be easily included in the equations of motion by using appropriate kinetic energy 
term in the Lagrangian , i.e., by using T instead of T rot and using p instead of p as a displacement 
vector. 

Using equations (59) and (70) the Lagrangian of the system is defined as 


L = T rot - V 


(71) 


For the purpose of convenience, L can be rewritten in the indicial notation as 


L = T r ot ~ V = \ 52 M HViPi ~ v (<l) 

The Euler-Lagrange equations for the system can then be derived from 


(72) 


— ( dL \ _ 2L 

dt \dpk) dp k 


(73) 


where, Fk are generalized forces from non-conservative force field. Evaluating the derivatives, 


and 




d_ (dL\ 
dt V<9p*J 


J2 M kjPj + 52 M kj pj 

3 3 


52 ***& + E 


dM ki . . 

nr** 


Also 


dL _ 1 v dMjj . . dV 
dpk 2 dpk PtPj dp k 


Thus the Euler-Lagrange equations can be written 


52 ^ + 521 


dMkj 

dp t 


2 d Vk >p ' p ‘ a Pk ~ Fl 

k = 1, ,...,n 


(74) 


(75) 


(76) 


(77) 


16 



By interchanging the order of summation and taking advantage of symmetry, it can be seen that 

v-~v r dMkj . 1 . f dM k j dM k { . . / r 7Q\ 

E { = o E { ~^r + -*r- ( 78 ) 

ij opi 2 dp, dpj 

Hence 

_ dM,. : 1 f)M:: — 1 ( dM kj dM kt 8M ijx . . ,_ 0 ^ 

iPiPj l. ,y j 


^ dM kj 1 dM, 3 . . t , v i rj. fc j 

Ei-^r 1 - = E + 


7j <9p, 2 dp* “2 dp, dpj dp k 


The terms 


^ _l f dM kj 8M ki dMij 

— 0 l a_ + J 


2 dp, dpj dp k 
are known as Christoffel symbols. Note that, for each fixed k, we have Cij k = C jk% . Also 

dV 


dp k 


— ^kjQj 


(80) 


(81) 


Finally, then Euler- Lagrange equations of motion can be written as 


^ ^ M kj Pj + ^ ^ C jjkPiPj T Dkjp I\.fzjPj — Ffc, ( k 1,2, rt) (82) 

3 iJ 

where D is the inherent structural damping matrix and Dp is the vector of nonconservative 
forces. 

In the equations (82), there are four types of terms. The first type of terms involve the 
second type derivative of the generalized coordinates. The second are quadratic terms in the 
first derivatives of p, where the coefficients may depend on p. These terms can be further 
classified into two types. Terms involving a product of the type p 2 are called centrifugal terms, 
while those involving a product of the type p t pj where i ^ j are called Coriolis terms. The third 
type are the ones which involve only the first derivative of the generalized coordinates and they 
are the dissipative forces due to the inherent damping. The fourth type of terms involve only p 
but not its derivatives. These arise from differentiating the potential energy. In the matrix-vector 
notation, the equations (82) are written in a compact form as 


M (p)p + C{p , p)p + Dp + Kp = F 


(83) 


The k y j - th element of the matrix C(p,p) is defined as 

n 

C kj = EcMPi 


1 = 1 


^ l ( dM kj< dM ki dMij 


= Ebi 


+ 


2 dp, dpj dp k 


)Pi 


(84) 
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Now, an important property of systems whose equations of motion are given by (83), is derived 


next. 

Theorem.- The matrix S{p,p) = M(p) - 2 C(p,p) is skew symmetric. 

Proof.- The kj - th component of the time derivative of the inertia matrix, M(p) is given by the 
chain rule as 


’ h 9p< 

Therefore, the fc/-th component of S = M — 2C is given by 


S kj = Affc, — 2Cki 


kj ~ t'-'kj 


= E 

t=i 

n 

= E 


dM kj _ 8M kj , dMx _ dM i} . 

dpi ' dp, dpj dp k \ 

dMij dM ki } 


ti [ 9 p„ dpj j " < 85 > 

Since the inertia matrix is symmetric, i.e., = M ju it follows from (85) by interchanging the 

indices k and j that 

£j k ~ &kj 


This completes the proof. 


3.6 Conclusions 

A generic mathematical model for a class of multibody flexible spacecraft was developed. A 
judicious choice of coordinate systems was made which allowed decoupling of translational and 
rotational dynamics of the spacecraft. This is a useful way of modeling spacecraft dynamics 
since, depending on the application, one can choose to use only rotational, or translational, 
or a complete rotational-plus-translational dynamic model of the spacecraft. Although only 
rotational model was presented, the Lagrangian formulation was done for a complete rotational- 
plus-translational dynamic model of the spacecraft. In the derivation of potential energy, an 
assumption was made that potential energy terms are only due to elastic deformations; however, 
as stated previously, the potential energy contributions from any other sources of deformations, 
such as thermal deformations, can also be included in the same way as elastic deformations. 
The model developed can be used for multibody spacecraft such as space-based manipulators, 
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multipayload platforms, satellites with flexible appendages, and many more. With minor modi- 
fications, the model can be used even for terrestrial robots. In summary, most spacecraft models 
can be obtained as special cases of the model developed in this paper. 
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Figure 1. Fixed reference and body axes systems 
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(b) Y - rotation 



(c) Z - rotation 


Figure 2. Basic rotations 
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